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Abstract 

Nonlinear spin motion in ferromagnets is considered with nonlinearity due to 
three factors: (i) the sample is prepared in a strongly nonequilibrium state, so that 
evolution equations cannot be linearized as would be admissible for spin motion not 
too far from equilibrium, (ii) the system considered consists of interacting electron 
and nuclear spins coupled with each other via hyperfine forces, and (hi) the sample 
is inserted into a coil of a resonant electric circuit producing a resonator feedback 
field. Due to these nonlinear ities, coherent motion of spins can develop, resulting 
in their ultrafast relaxation. A complete analysis of mechanisms triggering such a 
coherent motion is presented. This type of ultrafast coherent relaxation can be used 
for studying intrinsic properties of magnetic materials. 



PACS: 76.20. +q, 76.60.Es 
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1 Introduction 



There are several different types of spin dynamics in condensed matter, which can be 
distinguished according to whether the sample studied is in equilibrium, weak nonequi- 
librium, or strong nonequilibrium. Microscopic spin oscillations in equilibrium magnetic 
materials are related to magnons and are studied by scattering techniques, such as neutron 
[1] or light [2] scattering. Small deviations from equilibrium, caused by an alternating 
external field, are characteristic of resonance experiments, like electron spin resonance [3] 
or nuclear magnetic resonance [4]. However, when the initial state of a spin system is 
made strongly nonequilibrium, several types of spin relaxation can occur. If there are 
no transverse external fields acting on the spins, they relax to an equilibrium state by 
an exponential law with a longitudinal relaxation time TV When the motion of spins 
is triggered, at the initial time, by a transverse magnetic field, the relaxation is again 
exponential, but with a transverse relaxation time T 2 which is usually much shorter than 
TV 

A rather different relaxation regime from a strongly nonequilibrium initial state arises 
if the spin system is coupled to a resonator. This can be done by inserting the sample 
into a coil connected with a resonance electric circuit. Because of the action of the 
resonator feedback field, the motion of spins can become highly coherent resulting in 
their ultrafast relaxation during a characteristic time much shorter than T 2 [5]. This latter 
type of coherent spin relaxation from a strongly nonequilibrium state in the presence of 
coupling with a resonator is the most difficult to realize experimentally and to describe 
theoretically. Experimental difficulties have been overcome in a series of observations of 
this phenomenon for a system of nuclear spins in a paramagnetic matrix [6-10]. A theory 
of the coherent spin relaxation could be based on the phenomenological Bloch equations, 
but solely for the case when the process is triggered by a sufficiently strong coherent 
pulse thrust on spins at the initial time, so that spin interactions are of no importance 
and only the resonator field plays a role. However, the most interesting case is when the 
coherent relaxation develops in a self-organized way from an initially incoherent state, 
with no external coherent pulses triggering the process. For such a self- organized coherent 
relaxation spin interactions are of crucial importance. Then the Bloch equations become 
inapplicable and one has to resort to microscopic models. 

A microscopic approach for describing coherent processes in spin systems has been 
recently developed [11, 12] and applied to a system of nuclear spins interacting through 
dipole forces. It was shown that the main role in initiating self-organized coherent relax- 
ation is played by the anisotropic (so-called nonsecular) part of the dipole interactions. 

In the present paper we extend the microscopic theory of coherent spin relaxation 
[11, 12] to a much wider class of materials. We consider a rather general Hamiltonian 
including both nuclear as well as electron subsystems interacting with each other through 
hyperfine forces. The electrons can possess a long-range magnetic order as in ferromag- 
nets or ferrimagnets, and magnetocrystalline anisotropy is taken into account. A general 
investigation of strongly nonequilibrium nonlinear processes in realistic magnetic materi- 
als is of interest by itself and can also be useful for many applications. For instance, the 
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self-organized coherent relaxation, being quite different from other types of relaxations, 
may give additional information on intrinsic properties of magnetic materials. The ul- 
trafast relaxation can be employed for repolarizing solid-state targets used in scattering 
experiments [10, 13]. Coherent effects in spin systems, being similar to their coherent 
counterparts in optics [14, 15], could be used for analogous purposes but in another fre- 
quency region. For example, spin masers [16-18] can be realized. The sensitivity of the 
characteristic times of coherent relaxation to initial conditions could be used for creating 
ultrasensitive particle detectors [19]. 



2 Electron-Nuclear Spin Hamiltonian 

To make our consideration applicable to a wide class of magnetic materials, we take a 
rather general Hamiltonian 

H = H e + H n + Hi n t, (1) 

describing a realistic situation where the sample contains electrons with a Hamiltonian 
H e and nuclei with a Hamiltonian H n , their interaction being given by H int . The electron 
Hamiltonian is 

H e = — - ^2 JijSi • Sj — n e *y]B ■ Si, (2) 

i^j i 

where Jij is an exchange interaction; the indices i,j = 1, 2, N e here enumerate electrons; 

— * 

Si is a spin operator; /i e = g e fiB, with g e being the electronic gyromagnetic ratio and \ib, 

— * 

the Bohr magneton; B is a magnetic field. The nuclear Hamiltonian has the form [20] 
commonly accepted in the theory of nuclear magnetic resonance, 

4 = ^EE <%fl?lf -Hn^B-Ii, (3) 

ijtj a/3 i 

in which i,j = 1, 2, N n enumerate nuclei with dipole interactions 

< = #(<W-3<<-) (4) 

ij 

between each other, where fi n = g n f^N, Qn being the nuclear gyromagnetic ratio, fiN, the 
nuclear magneton, and r^- = \ fij \ , n^- = rij/rij, fij = — r}; the indices a and (3 label 
the components of Cartesian vectors (a, (3 = x,y,z); Ii is a nuclear spin operator. The 
general form of the hyperfine interactions [20, 21] between electrons and nuclei is 

H int = AJ2Si-Ti + lj2E ttfsfil (5) 

i iy^j a/3 

containing an isotropic contact part with an interaction intensity A and a dipole part 
with the interactions 

< = ^(^-3ng-ng). (6) 

ij 
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The total magnetic field is the sum 

B = H e z + H x e x , H x = H a + H, (7) 

of an external magnetic field in the z direction and of a transverse field including an 
effective field of a transverse magnetocrystalline anisotropy [22] and a feedback field H of 
a resonator. The longitudinal part of the magnetocrystalline anisotropy can be included 
into the external magnetic field H . 

In the preceding formulas we have used, for simplicity, the same indices, % and j, to 
enumerate electrons and nuclei, keeping in mind that for each particular case these indices 
run over different sets, so that for electrons i = 1, 2, N e and for nuclei i — 1, 2, N n . 
The corresponding electron and nuclear densities, p e = N e /V and p n = N n /V, where V 
is the volume of a sample, are, in general, different. 

The resonator coil is directed along the x-axis, so that the current induced in it is 
caused by the motion of the transverse magnetization 

M x = ±J2(p e <S?>+p n <I*>), (8) 

i 

where the angle brackets mean statistical averaging. The resonant electric circuit is char- 
acterized by a natural frequency u, ringing time 73, and quality factor Q, given by 

^im< ' K =m Q =t' (9) 

where L, C, and R are inductance, capacity, and resistance, respectively. The resonator 
feedback field is given [12] by the Kirchhoff equation 

^ + 2 l3 H + j H(r)dr = -4^^, (10) 


in which r] is a filling factor. 



3 Coupled System of Equations 

Spin dynamics is defined by the Heisenberg equations for the electron spin operators Sf 
and Sf, and for the nuclear spin operators if and if, where Sf = Sf±iSf, if = I?±ilf. 
These equations are coupled to each other by hyperfine interactions between electron and 
nuclear spins. Besides that, both types of spin motion are coupled with the resonator 
through the resonator feedback field, defined by equation (10), and the average transverse 
magnetization (8) expressed by means of the average spins. To describe the dynamics 
of spins coupled with each other as well as with a resonator, we need to derive the time 
evolution equations for the average electron and nuclear spins 

* = ^E<sr>, ^jfH<st>, (ii) 
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(12) 



The main steps of deriving these equations are the same as in Refs. [11, 12] except that 
now the Hamiltonian (1) is more complicated. 
We introduce the Zeeman frequencies 



and the anisotropy parameters 

where the Planck constant is set ft = 1. We use the notation 

£o = \ E (o<j < J j > < J / > +4 < J 7 >) > 



en < L > +2bu < If > +2cu < I- > 



in which 



= J {Aij ~ A^j - '-/.!';/) , Cij = - - /.Y'lj ) , 



A^f being given in Eq. (6). Also we write 



E [{aij - en) < I* > +Cij < If > +4 <I 3 > + 



. + 



V? = E [ 2 < J / > < J | >) + 



X - (etj < S~ > +2b tj < S + > +2cij < S] >) 



where 



6y = T(Cy -C^-2iC^ 



Qj = -{ct z - id; 



(13) 
(14) 

(15) 
(16) 



(17) 



(18) 



C°f being defined in Eq. (4). Equations (15) to (18) describe local random fields caused 
by spin fluctuations. In the uniform approximation, all these quantities would be zero, 
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because of the properties of dipole interactions. However, these local fields cannot be ne- 
glected, since they play a crucial role at the initial stage of spin relaxation. Therefore they 
must be retained and treated as local random variables. For other terms in the evolution 
equations, having a long-range nature in real space, one may employ the semiclassical ap- 
proximation. This approach of using for long-range terms the semiclassical approximation 
complemented by the stochastic quantization of short-range terms has been developed in 
Refs [11, 12]. 

Following these steps and taking into account the longitudinal, 71, Ti, and transverse, 
72, T 2 , attenuations for the electron and nuclei, respectively, we obtain the evolution 
equations for the electron spin variables (11) and nuclear spin variables (12). For the 
transverse and longitudinal electron spins we have, respectively, 

dx 

— = i(u e + - Co - As)x-i (a e + /i e H - £ - Au) z, (19) 
at 

dz % % 

— = -(a e + ii e H-i- Au) x* - - (a e + u e H - f * - Au*) x-^iz-a), (20) 
at 2 2 

where 71 and 72 are attenuation parameters and a is the stationary value of z. The 

variable x is complex, while z is real, so we should add either an equation for x* or for 

\x\. It is convenient to consider 



d \x\ 2 



—272 |x| 2 + i (a e + u e H — C ~ Au*) zx — i (a e + u e H — £ — Au) zx* . (21) 



dt 

In the case of nuclear spins, we obtain 
du 

— =i{oj n + iY 2 -<Po- Az) u-i(a n + u n H - ip - Ax) s (22) 

ds % % 

= 2 ( a ™ + ^ nH -V- Ax) u* — (a n + u n H -ip* - Ax*) u - T 1 (s - q) , (23) 

where Ti and T2 are the longitudinal and transverse attenuations, respectively, and q is 
the stationary value of s. In addition, we shall need 



d \uf 
dt 



= —2T 2 \u\ 2 + i (a n + u n H — tp* — Ax*) su — i (a n + fi n H — tp — Ax) su*. (24) 



All equations (19) to (24) contain the resonator feedback field H described by (10). 
The latter can be transformed [12] to the integral feedback equation 

t 

H = -47iT] J G(t - T)dM x (r), (25) 


expressed through a Stieltjes integral with the Green function 

73 

G(t) = (cosu 3 t sino^i) exp(— j 3 i) 

Us 
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and the differential measure dM x with M x defined in Eq. (8). Here the effective frequency 
is lj 3 = \]uj 2 - 73- 

The system of seven nonlinear equations (19) to (25) determines the dynamics of 
electron and nuclear spins coupled with each other as well as with a resonator. 



4 Scale Separation Approach 

Our aim here is to study the strongly nonequilibrium regimes of spin motion. This problem 
is different from considering the equilibrium properties of coupled electron and nuclear 
spins [23-26]. An additional complication, in our case, arises from the coupling of spins 
with a resonator by means of the feedback equation (25). To solve equations (19) to (25), 
we employ the scale separation approach [11, 12] which is a generalization of the averaging 
techniques of dynamical theory [27, 28] to statistical systems. 

To understand what different time scales exist for the system considered, we need 
to specify what small parameters we have. Since we have a sample coupled with a res- 
onator, some small parameters should appear by concretizing the corresponding resonance 
conditions, assuming the ringing width is much smaller than the natural frequency, 



^«1. 



UJ 



(26) 



The resonator natural frequency can be tuned either to the frequency of electron spin 
resonance uj e , so that 



UJ r , 



«1, 



uo — U) P 



(27) 



or to the frequency of nuclear magnetic resonance 

uj n = u n - Az eff , (28) 
in which z e ff is the longitudinal electron spin z averaged over the period 2n/u, so that 

A N 



«1, 



UJ 



U1 N . 



We assume that the external magnetic field H is sufficiently strong that 



«1, 



HeHeff 



«1, 



A 



«1, 



(29) 



(30) 



where H e ff is the resonator feedback field averaged over a period 2-k/uj. We also assume 



7i 



«1, 



72 



< 1. 



(31) 



Then from Eqs. (19) to (21) it follows that the variables z and |x| 2 are to be treated as 
slow compared to the fast variable x. Similarly, for nuclei we assume 



OLr. 



UJ N 



«1, 



UJ N 



«1, 



Ax, 



eff 



UJ N 



«1, 



(32) 
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where the subscript eff means again that the corresponding quantity is averaged over 
2ty/lj, and we keep in mind the usual inequalities 



«1, 



< 1. 



UN 



(33) 



Then Eqs. (22), (23), and (24) show that the variables s and \u\ are slow compared to 
the fast variable u. As the nuclear magnetic moment /i n is much smaller than that of an 
electron, /z e , we have the inequalities 



pn 



Ate 



<1, — <1, — <1. 

7i 72 



(34) 



Then, comparing Eqs. (19) and (22), we see that the variable u is slow compared to 
the fast variable x. And the comparison of Eqs. (21) and (24) tells us that |w| 2 is slow 
compared to the faster \x\ 2 . 

One more condition assumed is related to the local random fields (15) to (18). These lo- 
cal fields define the parameters of inhomogeneous broadening due to the electron-nuclear, 
T en , and nuclear- nuclear, T nn , interactions. These widths are assumed to satisfy 



U! N 



«1, 



U! N 



< 1. 



(35) 



Using conditions (30) and (32), we may simplify the feedback equation (25) to 



f „ dx „ du\ 



in which the parameters 



Pe 



2 VePe 
7T 7] 



2 pnPn 
71 T] 



(36) 



(37) 



characterize the effective coupling of the sample with the resonator. The details are given 
in the Appendix. Substituting into Eq.(36) the derivatives from (19) and (22), we find 

H = -2Re i(3 e [(u e - As - £ + ry 2 ) x + (Au + z] - 

- 2Re i@ n [(u n - Az - ip + iT 2 ) u + (Ax + ip) s] . (38) 

The feedback field (38) is to be substituted into equations (19) to (24). 

Then we solve Eq. (19) treating there all slow variables as quasi-integrals of motion. 
The solution reads 

x = (xq — x) exp | (iVt e — 72) £ j + x, (39) 

where x = x(0) and 



ijj P 



AS - ^0 - PePel2Z, 72 = 72 + VePe (^e ~ As - ^o) Z, 
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x=-^-{(a e -Au-C)z + ifi e p e [A (u* - u) + C - f] z 2 ) . 

After this, we solve Eq. (22), keeping the slow variables fixed and averaging the fast 
variables x(t) and z(t) to obtain x e ff and z e ff. Since x(t) is already known, we have 



x 



A e z/Q e , uj ~ c<j e 

0, CiJjv j 



(40) 



where Vt e = uo e — As — £ and A e = a e — Au — £. The solution of Eq. (22) is 

u = (u — u) exp | (iVt n — r 2 ) tj + m, 
where uq = -u(O) and 

O n = CJ„ - A (1 + /X„/3 e s) Z e // - ¥» - VfiPn^S, 
f 2 = T 2 + //„/?„ (w n -Az e ff- <p ) S, 
U = 7T ~ ^) s + »AW?e (C* - s ^e//] • 

The solutions (39) and (40) are to be substituted into the equations for the slow 
variables, with the right-hand sides of the latter equations being averaged over time and 
over random local fields according to the rule 



2Wo fit ^ ] 



dt 



with the stochastic measure m(ip) such that 

< <A) > = <</?> = < > = , 
< 6<aj > = < > = < > = « > = o , 
«£o» = «l£l 2 » = 7* 2 = rl, 



< y?o > = < M 2 > 



r 2 



r 2 +r 2 

nn 1 en 



The constants 7* and are the parameters of inhomogeneous broadening caused by 
hyperfine electron-nuclear dipole interactions and by nuclear dipole interactions. 
We introduce the effective coupling parameters 



2 PeHl^E ( Pn^nAs S 

g e = 7r r] I 1 + 



72 00 



_ 2 PnP n U N L 

g n = -K rj— I 1 



T 2 u 



PePeAz ef j 

p n p n uj N 



(41) 
(42) 



characterizing the strength of coupling between electron or nuclear spins, respectively, 
and the resonator. The effective frequencies are 



to E = u e - As, u N = u n - Az eff , (43) 

which are the frequencies of the electron spin resonance and nuclear magnetic resonance. 
Also, we define new slow variables for electrons, 

v = N 2 - ^ + A2 ' M|2 + 7 V , (44) 



and for nuclei 



where 



. |2 al + Yl + 5 2 2 /A _. 
w = \u\ - £ s, (45) 



N 



6 = V2n^^^z eff . (46) 
Accomplishing all these steps, we obtain from Eqs. (20) to (24) 

dz 

— = l2 g eV - ll ( z - a ), (47) 
^ = -2 l2 (l + g e z)v, (48) 

- = r 2 </ B «;-r 1 ( a -o, (49) 

^ = -2T 2 (l + g n s)w. (50) 

Equations (47) to (50) define the averaged motion of slow variables. We require their 
solutions, since all observable quantities can be expressed through them. 



5 Nuclear Spin Dynamics 

Equations (47) to (50) show that the electron spin dynamics is qualitatively similar to 
that of nuclear spins, but there are three main points distinguishing electron from nuclear 
relaxation. First, electron spin processes are usually much faster than nuclear processes 
[3], which is related to the fact that 72 3> T 2 . Second, as the electronic magnetic moment is 
three orders of magnitude larger than the nuclear magneton, electron spin motion is much 
less influenced by the presence of nuclei than the motion of nuclear spins by the existence 
of electrons. Third, the stronger influence of electrons on the motion of nuclear spins 
is caused by a long-range magnetic order that more readily occurs in electronic systems 
than in nuclear ones. Therefore nuclear spin dynamics is a little more complicated but at 
the same time richer than the dynamics of electron spins. 
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Suppose that the electron spins either were not perturbed at the initial time or, if 
perturbed, that fast electron processes have already been relaxed to their stationary state. 
Let us study the dynamics of nuclear spins that were initially prepared in a strongly 
nonequilibrium state. We denote the initial conditions for the nuclear spin variables as 



s(0) = s , w(0) = w , u(0) = u . 



(51) 



If relaxation of nuclear spins lasts for times of order Ti = T^ 1 , this would mean that 
coherent processes do not develop. Such a case would be of no interest for us, since our 
aim here is to investigate the fast coherent relaxation. Hence we shall consider times such 
that t <C rj" 1 . In this case, we may omit in Eq. (49) the term containing IY Then we 
can solve Eqs. (49) and (50) analytically obtaining for the longitudinal nuclear spin 



s = 



7o 

9^2 



tanh 



t-t 

T 



1 

9 



where, for the sake of simplicity, we write g = g n , and 

2 

I sech 2 



w 



7o 



t-t 



T 



Here 7 = r 1 is the relaxation width given by the equation 



7o 



'i + gso) 2 + g 2 w 



1/2 



with t being the relaxation time, and t the delay time, 



t n = — In 



7o - T 2 (1 + gs ) 



7o + r 2 (l + ^ ) 



(52) 



(53) 



(54) 



(55) 



According to Eq. (45), the modulus squared of the transverse nuclear spin is 

2 

I sech 2 



\u\ 



7o 



^r 2 



t - tr 



T 



, < + r l + s 2 2 

+ — t S. 



(56) 



A? 



If the coupling parameter g, is small, g <ti 1, then Eqs. (49) and (50) show that the 
relaxation of nuclear spins follows the standard exponential law with the relaxation times 
T\ = Ti 1 and T 2 = T^ 1 . This trivial regime is not interesting for us, so we concentrate 
attention on the case of strong coupling, when g » 1. Then the relaxation width (54) 
can be written as 

so 



7o = T 2 gJs^ + w + 



Sq + Wo, 



(57) 



For the delay time (55), we find 

9 (si + wo) + s - (1 + gs ) ^s% + w 



+ r ° i 
U = — In 



9 (si + w ) + s + (1 + gs ) y/sl + Wo 



(58) 
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For the relaxation time, after using again the inequality g >> 1, we have 



r = —7==. (59) 

A large value of the coupling constant g means, according to its definition in Eq. (42), that 
nuclear spins are strongly correlated with each other by means of an effective interaction 
through the resonator feedback field. As a result of this correlation they move coherently, 
which leads to the nonzero value of w \u\ 2 . Recall that, by definition (12), u = for 
incoherent spins. Coherent motion of the spins results in their ultrafast relaxation, which 
follows from Eq. (59) yielding To <C T2 when g ^> 1. This is why the characteristic time 
(59) can be called the coherent relaxation time. Notice also that the coupling parameter 
(42) is proportional to the number of nuclei N n , and so r ~ 1/g ~ l/N n . Such a 
dependence of the relaxation time, r ~ l/N n , on the number of radiators is typical for 
coherent processes that in optics are called superradiance [14, 15]. 

The relaxation characteristics, as is seen, essentially depend on the initial conditions 
(51). If Mo 7^ at the initial time, this implies that an initial coherence is imposed on the 
spins, which can be done by means of a short external pulse of a transverse field. When 
uq ~ 1, then we have wq ~ |«o| 2 , since the second term in Eq. (45) is small. In this case 
we get the regime of triggered relaxation [12]. 

A much more interesting question concerns how the relaxation of nuclear spins starts 
when no initial coherence is thrust upon the spins. This problem is also more important 
than consideration of the case when relaxation is triggered by external fields. For, when 
the relaxation is initiated not by external forces but by internal interactions, the relaxation 
characteristics, such as the delay time and the relaxation time, significantly depend on the 
parameters of these internal interactions. Therefore, the self-organized relaxation reflects 
(and provides information about) the intrinsic properties of matter. 

Let us consider the self- organized relaxation in the system of nuclear spins strongly 
coupled with a resonator. That is, we analyze the case corresponding to the conditions: 

m = 0, s 7^0, 0|s o |»l. (60) 

If uq = 0, then, according to relation (45), we have |wo| <C Sq- This permits us to simplify 
the relaxation width (57) getting 



7o = si 



g\s \ (2sq + wo) +e(2sl-woj] , (61) 



where e = sgn s . The delay time (58) becomes 

2 (1 - e) s 2 + w 



+ r ° 1 
to = y In 



and for the relaxation time (59), we get 



2(l+e) s 2 + w 



(62) 
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T o = — 7— 7* — — — 1 — r- (63) 
2g\s \ 3 T 2 g\s \ 

The delay time (62) strongly depends on the sign of the initial polarization of nuclear 
spins, e = sgn s . When this initial polarization is directed along the external magnetic 
field H , i.e., along the z-axis, then 



To 



to 



TO 



In 



w 



(e = l), 



and if the initial polarization is directed opposite to the external field, then 

t 



In 



4 s 2 
^*o 



w 



-1). 



(64) 



(65) 



In the case when the initial polarization is along the z-axis, Eq. (64) shows that 
to < 0, since \w \ <C s 2 ,. Then the function (53) quickly decreases starting from t — 0. 
This function being proportional to \u\ 2 , describes the degree of coherence in the 

motion of nuclear spins. In turn, the function \u\ is proportional to the power flow in the 
resonator circuit and, thus, is a directly measurable quantity [11, 12]. Hence, when w(t) 
quickly decreases starting from the initial time t — and \wo\ <C Sq, this means that no 
noticeable coherence develops in the system. 

In contrast, if the initial polarization of nuclear spins is directed opposite to the ex- 
ternal magnetic field, so that the system is prepared in a strongly nonequilibrium state, 
then, from Eq. (65), there is a positive solution for the delay time to > 0. In such a case, 
the function (53) increases from its initial value wq, reaching a maximum at t — to, when 



W (t ) = So, s (t ) = 



1 

9 



(66) 

- ^0; 



This means that a self-organized coherent pulse develops with a maximum at t 
which explains why t is called the delay time. 

Combining both the cases, (64) and (65), into one and substituting the expressions 
for T and w , we obtain 



^9 so 



In 



a' 



+ r 2 +r 2 +5 2 

1 x en 1 nn 1 " 



4w 2 



N 



(67) 



This is the central formula for analyzing what internal microscopic mechanisms are re- 
sponsible for the self-organized development of coherent relaxation. Each of these internal 
mechanisms is related to the corresponding parameter entering formula (67). Among such 
internal causes that may trigger self-organized coherence, we have the transverse magne- 
tocrystalline anisotropy, a n ; the dipole part of the hyperfine interactions, characterized 
by the parameter T en ; the dipole interactions between nuclear spins, T nn ; and the param- 
eter 5 defined in Eq. (46), which is due to the simultaneous existence of the hyperfine 
interactions, of coupling with a resonator, and of the magnetization of electron spins. 
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The relaxation time (54), for the case of self-organized relaxation, when u = 0, reads 



r = T 2 



2 , 2 2 "n 



2 i r2 i r2 i X21" 1 / 2 



(l+gs y + g 2 s 2 



a 2 +T 2 + r 2 +5 



(68) 



This demonstrates that the value of the relaxation time (68) depends mainly on the 
strength, g, of coupling with a resonator. In this way, the delay time (67) and the 
relaxation time (68) are related to different characteristics of the system considered. 

In order to decide what kind of interactions, direct dipole interactions between nuclei 
or hyperfme interactions, influences more the values of the characteristic times, we should 
compare the corresponding widths T nn ~ p n p n an d Y en ~ pp e Pn, where p ~ min {p e , p n }. 
There are two limiting cases. The first is when p n < p e , and then 

L """ / '" 10 ^ (Pn<Pe)- (69) 



Ten p 



e 



Hence, when the density of nuclei is lower or comparable with that of electrons, nuclear 
dipole interactions are negligibly small compared to the hyperfme interactions between 
nuclei and electrons. Another case is when the density of nuclei is much higher than that 
of electrons; then 

Tnn PnPn Pn n — 3 



^10"' (p„» Pe ). (70) 

J- en Pepe Pe 

Thus, the nuclear dipole interactions become stronger than the hyperfme interactions only 
when the density of nuclei surpasses by three orders of magnitude the density of electrons. 

In equilibrium theory, the influence of hyperfme interactions is often modelled by the 
effective Suhl-Nakamura forces directly acting between nuclear spins [21]. These forces 
are responsible for the appearance of nuclear spin waves corresponding to well defined 
excitations, even at those temperatures where the nuclear spins are completely disor- 
dered. The underlying cause of the formation of the nuclear spin waves is the existence of 
magnetic long-range order in the electronic subsystem, which defines both the long-range 
interaction radius of the Suhl-Nakamura force and its existence as such. The effective 
Suhl-Nakamura force describes an indirect interaction of nuclear spins through magneti- 
cally ordered electrons [21]. 

For strongly nonequilibrium processes, such as those considered in this paper, the role 
of the magnetic order of the electrons is essentially different. This order does strongly 
influence several important characteristics. For instance, in addition to the usual shift of 
the nuclear magnetic resonance frequency cun = u n — Az e ff, it leads to the appearance 
of parameter (46) playing for nuclei the role of an additional inhomogeneous width. Nev- 
ertheless, even if this order is absent, so that z e ff — * 0, coherent nuclear spin relaxation 
can exist. When z e ff — > 0, the Suhl-Nakamura force is not well defined, but the hyper- 
fine interactions do not stop existing. These interactions define the width T en , which is 
not zero even if the electronic magnetic order is absent. Thus, the presence of hyperfme 
interactions is already important, even when there is no long-range magnetic order, when 
the Suhl-Nakamura force and nuclear spin waves are not well defined. 
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However, it is worth emphasizing that the appearance of electronic magnetic order 
can strongly change the values of the characteristic parameters. Thus, one of the most 
important parameters is the effective coupling (42) describing the coupling of nuclear 
spins with a resonator. The value of this parameter essentially depends on whether z e ff 
is zero or not. The appearance of magnetic order in the electronic system can change 
the value of the parameter (42) by several orders of magnitude which, in turn, drastically 
changes the values of the delay time to and the relaxation time r . 

6 Characteristic Numerical Values 

To understand better the role of different factors in the coherent relaxation of nuclear spins 
and the magnitudes of the related characteristic parameters, let us now make numerical 
estimates. We take the values of parameters that are typical of many ferromagnetic 
materials in which one usually studies nuclear magnetic resonance and nuclear spin echo 
[21, 29, 30]. These can be pure materials, such as Co, or various ferromagnetic alloys 
and compounds [29, 30]. Since ferrimagnets are often treated by effective ferromagnetic 
models, ferrimagnetic materials, such as MnFe203, are also included here [31]. 

For the characteristic magnetic fields and the corresponding frequencies we have the 
following values. The contact hyperfine field Ha = A/p n ~ 10 5 G, the related frequency 
uja = A/h ~ 10 9 s -1 . The hyperfine field is smaller than the electron exchange field 
Hj = J//i e ~ 10 6 G, the corresponding frequency being uj = p e Hj/h ~ 10 13 s _1 . 
However, both these fields are important for nonequilibrium processes in the nuclear spin 
system, although the role of these fields is different. The hyperfine field acts directly on the 
nuclear spins, and the exchange field influences nuclear spin relaxation indirectly, through 
the formation of magnetic order in the electronic subsystem. If we take an external 
magnetic field H ~ 10 4 G, then the Zeeman frequencies (13) are u e = p e H /h ~ 10 11 
s _1 and uo n = p: n H /h ~ 10 8 s _1 . The magnetocrystalline anisotropy field H a < 10 3 
G, depending on the particular structure of matter. The anisotropy parameters (14) are 
a e = p e H a /h < 10 10 s _1 and a n = p n H a /h < 10 7 s -1 , respectively. The longitudinal 
widths 71 and I\ can vary within rather wide intervals, but usually 71 <C 72 and Ti <C IV 
For the transverse widths we may take the estimates 72 ~ p e p? e /Ti and T 2 ~ p n p\/Ti. This, 
with p e ~ 10~ 20 erg/G, /i n ~ 10~ 23 erg/G, and p e ~ p n ~ 10 23 cm -3 , gives 72 ~ 10 10 s _1 
and T 2 ~ 10 4 s -1 . The value for 72 is to be treated as the upper limit, since the density of 
electrons is usually less than 10 23 cm -3 , being, for instance, 10 22 cm -3 for typical ordinary 
metals [32]. In the case when the considered electrons are related to impurity ions inside 
an insulator, as in Refs [6-10], then their density can be p e ~ 10 20 cm -3 , resulting in 
72 ~ 10 7 s~ 1 . The estimated value of T 2 is in agreement with experimental measurements 
[30]. For the resonator ringing width, we may take a typical experimental value of 73 ~ 10 6 
s -1 . Then 73/0^ ~ 10~ 5 and 73/u^ ~ 10~ 2 , so that inequality (26) is satisfied. 

Since 71 72 and 72/we < 10 _1 , condition (31) is valid. The nuclear magnetic 
resonance frequency (28) is ~ u n ~ 10 8 s™ 1 if z e ff = 0, that is if the magnetic order 
is absent, and if z e ff 7^ 0, then un ~ 10 9 s -1 . Hence, r 2 /u;7v ~ r 2 /cj n ~ 10~ 4 , when 
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z e ff = 0, and T 2 /ujn ~ 10~ 5 for a ferromagnetic material with z e ff 7^ 0. This, together 
with Y 1 <C r 2 , shows that condition (33) holds true. 

For electrons, a e /u e < 10~ 4 and uj^/uj e ~ 10~ 2 , while for nuclei, a n /u^ < 10~ 2 
and x e ff = 0. The resonator feedback field (38) is H e ff ~ P e ^A- For electrons, with a 
resonator natural frequency close to the electron spin resonance frequency, cu ~ u e , we 
have fi e H e ff/u e ~ 10~ 3 , and for nuclei, when u ~ un, we find p n iY e ///^iv ~ 10~ 2 . Thus, 
all inequalities in equations (30) and (32) are valid. Since \p n /Pe\ ~ 10~ 3 and Ti/7i and 
are of the order of p n /i 2 /p e/ u 2 ~ (p n /p e ) 10~ 6 , the inequalities in Eq. (34) hold 
true if p n and p e are not drastically different. Conditions (35) are also satisfied, since 
T en /u N ~ p n ^ 2 /p e ^N ~ 10~ 2 and T nn /uj N ~ T 2 /uj n ~ 10" 5 . 

Among the parameters defining the characteristic times of the coherent nuclear spin re- 
laxation, we have the anisotropy parameter a n < 10 7 s _1 , the inhomogeneous broadening 
due to hyperfine dipole interactions, T en ~ pp e pn with p = min {p e ,p n }, the inhomo- 
geneous broadening due to nuclear dipole interactions, Y nn ~ p n p 2 ~ 10 4 s _1 , and the 
parameter 5 is given by Eq. (46), from which 5 2 ~ \§~ 2 Y 2 en z 2 e ^ . The width T en , ac- 
cording to Eqs. (69) and (70), is always larger than T nn , provided that the density of 
electrons p e is not three orders smaller than the density of nuclei p n , which follows from 
the relation T en ~ 10 3 (p e /p n ) r n „. For example, if p e ~ p n , then T en ~ 10 3 T nn . If we 
take p e ~ 2 x 10 20 cm -3 and p n ~ (5 x 10 22 — 10 23 ) cm -3 , as in the experiments [6-10], 
then T en ~ (1 — 10)r nn . In this way, the hyperfine width T en is usually larger than T nn 
and always larger than 5, although Y en may be comparable with a n , if p e p n . When 
Pe ~ Pn, the largest parameters among those considered above are a n and T en ~ 10 7 
s -1 . In such a case, other parameters entering additively with these can be omitted. For 
instance, the delay time (67) may be written as 



An important parameter entering the expressions for the characteristic times and 
influencing the behavior of solutions is the coupling parameter g, which is drastically 
different for the case when there is magnetic order in the electron system compared to 
the case when the magnetization is absent. For the latter case, when z e ff = 0, we have 
g ~ 10. When a ferromagnetic material is considered, so that z e ff ~ 1, then the second 
term in Eq. (42) can become much larger than the first one. Thus, for p e ~ p n , we have 



Therefore, the coupling parameter g can be increased by three orders by the presence 
of electron magnetization, reaching the value g ~ 10 4 . The system of magnetized electrons 
acts as an additional resonator strongly strengthening the coupling between the resonance 
electric circuit and nuclear spins. 

To evaluate the characteristic values of the delay time (67) and relaxation time (68), we 
made calculations for several ferromagnetic materials with typical parameters taken from 



In 




en 



p e p e Az e ff p e 

rs_/ 

PnPn^N Pn 



10 3 . 
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Ref. [30]. In our formulas we take the filling factor rj — 1, consider the purely resonance 
case, when u = \oon\, also take z e ff — |, p e — p n and assume that the transverse 
anisotropy parameter is small compared to T en . We analyze the case of purely self- 
organized coherent relaxation when at the initial time nuclear spins are polarized against 
the external magnetic field, so that s = —I, where I is an absolute value of a nuclear 
spin, and there is no initial coherence imposed upon the system, so that u = 0. The high 
initial polarization of nuclear spins can be achieved by the dynamic nuclear polarization 
technique. The transverse relaxation time T2 = T^ 1 can be measured by several methods, 
e.g., the two-pulse echo technique or the single-pulse echo technique [30], of which the 
former is likely to be more reliable. Our results are presented in Table I. 

7 Discussion 

Nonlinear spin dynamics is considered for ferromagnets consisting of electron and nu- 
clear subsystems coupled through hyperfine forces. The sample is prepared in a strongly 
nonequilibrium initial state. In addition, the ferromagnetic sample is considered inserted 
into a coil of a resonant electric circuit. All this makes the spin dynamics highly nonlin- 
ear. The evolution of the system is described by seven nonlinear equations, six of which 
are differential equations for electron and nuclear spins and one equation is an integro- 
differential equation for the feedback field of the resonator. These are solved by using the 
scale separation approach [11, 12]. It is shown that due to the resonator feedback field 
an ultrafast coherent relaxation of spins can occur. The system of magnetized electrons 
serves as an additional resonator for the nuclear spins, significantly enhancing the effective 
coupling of the nuclear spins with the resonator circuit. Such an enhancement can reach 
three orders of magnitude, as compared to the coupling in a paramagnetic material. 

The ultrafast coherent relaxation of nuclear spins may be either triggered by an ini- 
tial pulse or can be self-organized. The latter case is the more interesting, since then all 
relaxation characteristics, such as the delay time and relaxation time, depend strongly on 
the values of the internal parameters. The most important such parameters, starting the 
process of self-organized coherent relaxation and, therefore, defining the main relaxation 
characteristics, are the transverse magnetocrystalline anisotropy and the dipole hyperfine 
interactions. If the density of electrons is more than three orders of magnitude lower 
than the density of nuclei, then the direct nuclear dipole interactions also become impor- 
tant. An interesting extension of the present approach could be the inclusion of external 
alternating magnetic fields, as has been done for nuclear magnets [33-35]. 

By studying the peculiarities of the coherent spin relaxation, it is possible to extract 
information on the intrinsic properties of magnetic materials. This especially concerns 
the regime of self-organized coherent relaxation, whose characteristics are very sensitive 
to the values of the parameters of the material studied. By observing a coherent pulse 
in the power flow, one can measure, with a very high precision, the delay time to. The 
latter can be accurately measured because it exactly corresponds to the maximum of the 
function w(t), which is proportional to the power flow in the resonant circuit [12]. If the 
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delay time to is measured experimentally, then, inverting (67), one may find the sum of 
a 2 n and 1^ = T 2 en + Y 2 nn as 

«n + r * = 4uj 2 n exp (2gs ^ 

where the inequality 5 <C T en is taken into account. This relation, when a n is known from 
other experiments, makes it possible to define the inhomogeneous broadening T*. As is 
mentioned in Sec. VI, one usually has a n < T*. The transverse-anisotropy parameter a n 
depends on the orientation of the sample with respect to the external magnetic field. It 
is possible to choose an orientation such that a„ < r». Then we obtain a simple formula 
giving the inhomogeneous broadening 



I\ = 2un exp (gs 



T 2 



in terms of the known values of the nuclear magnetic resonance frequency un, the coupling 
parameter g = g n in Eq. (42), the initial nuclear polarization sq, the transverse relaxation 
time Tii an d the measured delay time to- It follows from the analysis of Sec. V that 
t > requires s Q < 0. Thus the value s t in the formula for is negative, making 
T* <C uj n . The exponential dependence of on the delay time t makes the value of 
very sensitive to to- 

Another possibility for exploiting the effect of coherent spin relaxation is its sensitivity 
to initial conditions, in particular, to the initial amplitude of the transverse spin, \uq\. The 
latter, in order to influence the delay time (58) and relaxation time (59), should be such 
that 

2 > o£ + T2 g2 1Q _ 4 

This is always small, and can be made arbitrarily smaller by reducing sq. Hence, we 
conclude that even quite weak external pulses, resulting in nonzero \uo\, can trigger the 
process of coherent relaxation. For example, from Eqs. (59) and (45) we get 

|«o| 2 = 

This allows us, by measuring the coherent relaxation time tq, to find the initial amplitude 
\u \. The sensitivity of coherent spin relaxation to initial conditions could be employed 
for creating ultrasensitive detectors of weak external pulses. In turn, this can be used to 
construct sensitive particle detectors [19]. 

In conclusion, the main results obtained in this paper can be summarized as: 

(i) A theory of nonlinear spin dynamics is developed for the systems of electron and 
nuclear spins coupled with each other through hyperfine forces and also coupled to a res- 
onator electric circuit. This essentially generalizes the previous consideration of nonequi- 
librium nuclear magnets [11-13] to a much wider class of materials, having long-range 
magnetic order. 
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(ii) The very complicated set of nonlinear differential equations is solved by invoking 
the scale separation approach [11, 12]. It is important that, because of the existence of 
small parameters resulting in different time scales, the motion of electron and nuclear 
spins can be effectively separated, as is seen in Eqs. (47)-(50). 

(iii) The effect of a strong coupling-parameter enhancement, due to the presence of 
magnetic order in the electronic subsystem, is described. The effective nuclear coupling 
parameter can be enhanced by three orders of magnitude, which makes relaxation re- 
ally ultrafast, with the relaxation time (63) becoming smaller than T 2 by four orders of 
magnitude. 

(iv) The nature of all main intrinsic mechanisms triggering the self-organized coherent 
relaxation, that occurs in the absence of external pulses, is elucidated. These mechanisms, 
defining the delay time (67), are the electron-nuclear interactions through hyperfine dipole 
forces, nuclear dipole interactions, and magnetocrystalline anisotropy fields. 

(v) Two types of applications are discussed. One type concerns the investigation of 
the internal properties of the materials by measuring the delay time (58) and coherent 
relaxation time (59). Another type utilizes the sensitivity of these characteristic times to 
the initial conditions, giving the possibility of employing coherent spin relaxation for the 
ultrasensitive detection of weak external pulses. 
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Appendix. Feedback Field 

According to the scale separation approach [11, 12], we first consider Eq. (19) for the 
fastest variable x, treating there the slow variables s and z as quasi-integrals of motion. 
Because of the second inequality in Eq. (32), we can, in a first approximation, omit the 
term containing H, whereupon the solution of (19) is 

with Q e = u e — As — £ and A e = a e — Au — £. In Eq. (22) we keep, as a quasi-integral 
of motion, the slow variable s. Because of the second inequality in Eq. (32), we again, to 
the first approximation, may omit the term with H. Averaging, in the right-hand side of 
Eq. (22), the electron variables over the period 2n/u, we get an approximate equation 

du . ,„ .„ x . . 
— ~ % [il n + iY 2 ) u - iA n s, 
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where Q n = u n — Az e ff — y?o and A n = a n — Ax e ff — (p. We find 



A n s 



u ~ \u — 



exp {i (Q n + iT 2 ) t} + 



A n s 



Substituting the approximate expressions into Eq. (25), we come to the form 

H = -2Re 



_ , s doc , , du 



in which 



A (0 = ^ r f]~~~~ [exp (* e f) - 1) , A. W = [^P (S n t) - 1] , 



5 e = i (u - iu e + As + £ ) + 72 - 73, 5 n = * - w« + Az + <p ) + T 2 - 73. 

The functions /3 e (t) and /3 n (t) do not vary much during the period 2n/u, so we may 
replace them by their averages over this time, 

2tt/lu 







2n5 e 
3a; 



^ 



2ix5 n 
3uj 



Omitting here the small terms 5 e /uj and 5 n /u>, we obtain Eq. (36). 
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Sample 


Nucleus 


/ 




T 2 


To 


to 








(10 9 Hz) 


(10~ 4 s) 


(10~ 8 s) 


(IO" 8 s) 


Lio.5Fe2.5O4 


57p e 


1/2 


0.47 


40 


88.2 


400 


Mno.5xSbo.49O4 


123 Sb 


7/2 


1.31 


1.70 


0.54 


2.98 




55 Mn 


5/2 


1.61 


0.60 


0.26 


1.53 


NiMnSb 


55 Mn 


5/2 


1.88 


0.95 


0.42 


2.49 


NiMnSi 


55 Mn 


5/2 


2.01 


0.60 


0.26 


1.59 


Co 2 MnSi 


59 Co 


7/2 


0.91 


0.38 


0.12 


0.62 




55 Mn 


5/2 


1.59 


0.80 


0.35 


2.03 


Co (fee) 


59 Co 


7/2 


1.37 


0.30 


0.09 


0.53 


Co (hep) 


59 Co 


7/2 


1.38 


0.65 


0.21 


1.15 



Table 1. The characteristic parameters related to the self-organized coherent nuclear 
spin relaxation in several ferromagnetic materials. 
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